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■ We discuss in this paper the subleading contact interactions, or counterterms, in the Po 

channel of nuclcon-nucleon scattering up to 0(Q 3 ), where, already at leading order, Wein- 
berg's original power counting (WPC) scheme fails to fulfill renormalization group invariancc 
due to the singular attraction of one-pion exchange. Treating the subleading interactions 
as perturbations and using renormalization group invariance as the criterion, we investi- 
gate whether WPC, although missing the leading order, could prescribe correct subleading 
^ _ counterterms. We find that the answer is negative and, instead, that the structure of coun- 



terterms agrees with a modified version of naive dimensional analysis. Using Pq as an 



■ example, we also study the cutoffs where the subleading potential can be iterated together 
> ' 

■ with the leading one. 
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Since its proposal [1 



Since its proposal |l[ and first implementation [2J, the effective field theory (EFT) description of 



00 

nuclear forces has enjoyed significant phenomenological successes 



Based on naive dimensional 



analysis (NDA), Weinberg's original power counting (WPC) scheme has one crucial assumption: 
each derivative on the Lagrangian terms is always suppressed by the underlying scale of chiral 
EFT, Mhi, where My ~ m a with m a the mass of the a meson. Due to the lack of bound states 
or resonances near threshold (with the delta-isobar integrated out), this assumption turns out 
to be quite applicable in the single-nucleon sector. However, the nonperturbative nature of few- 
nucleon systems could make infrared mass scale {M\ Q ~ 100 MeV) — such as the pion decay constant 
/tt ~ 92 MeV and the pion mass m n ~ 140 MeV — to enhance NN contact interactions relative to 
WPC. In fact many works have pointed out that WPC is inconsistent from the point of view of 



renormalization group (RG) invariance 



111 ], Here renormalization is associated with ultraviolet 



(UV) momenta in the Lippmann-Schwinger (LS) equation rather than with the UV divergence of 
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perturbative, NN irreducible diagrams. 

One mechanism to upset WPC is the singular (diverging at least as fast as 1/r 2 ) attraction of 
the tensor force of one-pion exchange (OPE), —1/r 3 at r — > 0, in certain S, P, and D channels 
where OPE needs full iteration: 3 Si - 3 Di, 3 P , 3 P 2 - 3 F 2 , and 3 D 2 . The singular attraction of 
OPE requires a counterterm at leading order (LO) in these channels. However, the P or D wave 
counterterms are counted as subleading in WPC since they have at least two derivatives [8(. In 
other words, the LO counterterm of, say, 3 Po, is enhanced by 0(M^/M?). We denote different 
orders of the EFT expansion by its relative correction to the LO. Thus, the next-to-leading order 
(NLO) is labeled by 0{Q/M hi ) or 0(Q) for short, and next-to-next-to-leading order (NNLO) by 
0(Q 2 /^hi) or 0(Q 2 ), and so on. 

The question remains as to how to modify WPC at subleading orders. The model study in 
Ref. fl^ suggests, referred to in the paper as modified NDA (NDA), that in the case of the 
LO long-range potential being singular and attractive, the subleading counterterms are enhanced 
relative to NDA by the same amount as the LO counterterms. As a first study, we use 3 Po in the 

? T. Note that different approaches 



paper to investigate whether NDA is applicable in nuclear E 
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and that perturbative pion 



toward renormalization of nuclear forces are offered in Refs. 
theory leads to a very different renormalization program [3- 

The resummation of ladder diagrams in lower partial waves can be achieved by solving the LS 
equation, with the schematic form, 

T = V + J VGT, (1) 

where the potential V consists of short-range counterterms Vs and long-range pion-exchange Vl, G 
is the Schrodinger propagator and A is a (sharp) momentum cutoff. Being singular and attractive 
in 3 Po, OPE demands promotion of the leading 3 Po counterterm, 

(»p„|vi»)|'p ) = a PllJ ,'p~^, (2) 

m N M lo 

an 0(M^/M 2 ) enhancement over WPC [8(, where ^/m^ is a common factor associated with 
nucleonic loops. For a comparison, the leading S-wave counterterms scale generically as Cs ~ 
4Tr/(m N M lo ). 

It will prove useful to have the short-range behavior of the LO 3 Po wave function. For kr <C 
1 < Ar where k is the center-of-mass momentum, the LO wave function ■0l°^( r ) can be solved for 



in powers of k l 



20, 



2l|, up to a normalization factor, 



4 0) (0~ 



AV 
r 



u + k 2 r\l- Ul + O(k A ) 



(3) 
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where A = ^"^f with = 1-29 the nucleoli axial charge, uq and U\ are oscillatory functions 
in terms of r/X and <j) with amplitudes ~ 1, where (j) is the phase between the two independent 
solutions and is related to C3p Q . 

Before two-pion exchanges (TPEs) are accounted for, which provide 0{Q 2 ) or higher corrections 
to OPE, there might exist a nonvanishing 0(Q) manufactured by inserting subleading contact 
interactions alone. In fact, Ref. [9J rated the triplet P-wave subleading contact interactions as 
0(Q 3 / 2 ), a lower-order contribution than TPE. 

We argue that a nonvanishing 0(Q) does not seem to be necessary. The study of Ref. 21] 



suggests that the residual cutoff dependence of the renormalized LO 3 Po amplitude is C(A~ 5//2 ). 
Had this residual cutoff dependence vanished slower than 0(A~ 2 ), our ignorance of short-range 
physics would have been larger than 0(Q 2 /M^ ; ), that is, larger than what TPE could compensate. 
If this happened, one would have had to consider an 0{Q /M^) correction to the LO amplitude 
induced by inserting a four-derivative 3 Po counterterm, before accounting for TPE. However, the 
rather small residual cutoff dependence of the LO 3 Po amplitude does not ask for a nonvanishing 
0(Q). 

With O(Q) vanishing, the 0(Q 2 ) and 0(Q 3 ) amplitudes, T^> and consist, respectively, 
in one insertion of 0(Q 2 ) or 0(Q 3 ) potentials, 

T (2,3) = y(2,3) + y A y(2,3) GT (0) + J r (0) Qy (2, 3) 

+ / / T^GV^GT^. 



Since this is equivalent to first-order distorted wave expansion, one can evaluate the "superficial" 
divergence of one insertion of (leading TPE) and V? (subleading TPE) before any coun- 

(2 3) 

terterm is considered, by investigating the short-distance behavior of the matrix element of V£ 
between the LO wave functions. This is facilitated by the short-distance behavior of the LO wave 
function © and TPEs, V^ 2) ~ 1/r 5 and ~ 1/r 6 . With a radial coordinate cutoff TZ ~ 1/A, 



we arrive at 



= (^(°)|yf V (0) > ~ j drr 2 \^°\r)\ 2 ^ 

~ a (A)A 5/2 + /3 Q (A)k 2 + C(fc 4 A" 5 / 2 ), (5) 

T (3) = <^°>|V?V (0) > ~ / drr 2 \^°\r)\ 2 \ 

J~l/A r 

~ qi(A)A 7/2 + /3i(A)A/c 2 + 0(k A A~ 3/2 ) , (6) 
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where ao,i(A) and /3o,i(A) are oscillatory functions diverging slower than A. Their exact forms can 



be evaluated [22j but are not crucial for our discussion. 

It is not necessarily true that one must use two counterterms to subtract the two divergent 
terms in Eq. ([5]) or ([6]). In fact, WPC prescribes only one counterterm up to 0(Q 3 ), and with the 
nonperturbative treatment it did provide a good fit to partial wave analysis (PWA) for a moderate 



range of cutoffs 



a. 



3H5J]. We first consider the counterterms prescribed by WPC, 



( 3 Po|# 3) | 3 Po) = cg o 3) ^, (7) 



where we split C3p into three pieces with C^p o determined at LO. This splitting reflects the fact 
that the value of the "bare" C3p could be modified at each order by, e.g., the short-range core of 
TPE, but the number of physical, short-range inputs is still one. 

The other scenario is to provide an equal number of counterterms as the divergent terms in 
Eqs. © and j6]): 

( 3 Po|^ 2 ' 3) | 3 P ) = C 3 ( P ' 3) p'p + Aip^ p'p(p' 2 + P 2 ) ■ (8) 

Thus, up to C(Q 3 ), every 3 Po counterterm gets enhanced relative to WPC by the same amount, 
0(M£/M£). This is exactly what to be expected based on NDA Q. 

There are a few versions of TPEs in the literature with slight differences in how double counting 
b »o id e d flfl Q. For define W e USe t he version b Re, fl. de^ess TPE e XpreS8i o„ S 
with dimensional regularization. We adopt the following low-energy constants for the v = 1 ttttNN 
seagull couplings (GeV -1 ): c\ = —0.81, C3 = —4.7 and C4 = 3.4 [3J. 

We first test the prescription given by WPC (|7|). At each order, Cs\> is determined such that 



the phase shift at = 50 MeV agrees with the Nijmegen PWA [24]] . Note that extra cares are 
needed to convert the subleading T-matrix into phase shifts, as shown in the Appendix. 

Figure [U shows the resulting phase shifts as functions of A at 7] a b = 40, 80, and 130 MeV. 

At 0(Q 2 ), the oscillatory cutoff dependence becomes more evident as the energy increases. This 

(2) 

is consistent with the superficial divergence ([5]); while the «o term is taken care of by the 
amplitude of the oscillation of (3o(A)k 2 is left intact and grows as the energy increases. Moving 
on to C(Q 3 ) we find more drastic, oscillatory cutoff dependence with a visibly growing amplitude. 
This is due to the factor of A that accompanies the oscillatory /?i(A) in Eq. ([6|). 



With the NDA counting ([8]), we need two physical inputs to determine the values of C3p and 
Z?3p . We fit them to reproduce the PWA at Ti a b = 20 and 50 MeV. The plateaus in Fig. [2] clearly 
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FIG. 1: (Color online) With the subleading counterterms ©, the 0(1), 0(Q 2 ) (upper row) and 0(Q 3 ) 
(lower row) 3 Pq phase shifts as functions of the momentum cutoff at Ti a b = 40 (a), 80 (b), and 130 (c) McV. 
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FIG. 2: (Color online) With the subleading counterterms ©, the 0(1), 0{Q 2 ), and 0(Q 3 ) 3 P EFT phase 
shifts as functions of the momentum cutoff at Ti a b = 40 (a), 80 (b), and 130 (c) McV. 

show the RG invariance of the power counting (JSj) , where the phase shifts are plotted as functions 
of A at given Ti^. 

In Fig. [3l the EFT phase shifts are plotted as function of energy. The fit is refined by employing 
more PWA points (?i a h =25, 50, 75, and 100 MeV) in the fitting procedure. We see that both 
0(Q 2 ) and C(Q 3 ) are in good agreement with the PWA. 

Although the perturbative treatment of WPC does not lead to cutoff independent results, the 
nonperturbative treatment does seem to fulfill RG invariance [l^]. It is therefore instructive to 
compare the following three scenarios for 3 Po: (i) the perturbative (Pert-CD) and (ii) nonpertur- 
bative calculations with the modified power counting ([8]) (Iter-CD), and (iii) the nonperturbative 
calculation with WPC (Iter- WPC). 
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FIG. 3: (Color online) With the subleading counterterms (j8]) and the improved fitting procedure, the 3 Po 
EFT phase shifts at <D(Q 2 ) and 0(Q 3 ) as function of 7I ab for A = 1500 MeV. 



- - Iter-CD 
Pert-CD 
■- Iter-WPC 
■ ■ PWA 




50 100 150 200 50 100 150 200 



T lab (MeV) 



T |ab (MeV) 



FIG. 4: (Color online) The 3 Po phase shifts by Iter-CD, Pert-CD and Iter-WPC as functions of laboratory 
energy, (a) is plotted with A = 400 MeV and (b) with A = 1200 MeV. 



Shown in Fig. U] are the 3 Po phase shifts calculated at 0(Q S ) with the aforementioned three 
schemes, where the fit of C3p and Dsp in Eq. ([8]) is performed with the PWA inputs up to 
7|ab = 50 MeV. At the lower end of cutoffs (as exemplified by A = 400 MeV), the three curves 
differ drastically from each other above Xi a b = 50 MeV, with the Pert-CD curve agreeing somewhat 
better with the PWA. As the cutoff goes higher (exemplified by A = 1200 MeV), the difference 
between Iter-CD and Iter-WPC becomes smaller and eventually vanishes, in accordance with the 



finding of Ref. 25(. We note that this does not necessarily mean that the fitting drives Da Po to 0. 
Rather, the quality of the fit is not sensitive to Z?3p A 2 /C3p when the ratio is tuned from to 1. 
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In summary, we conclude: 

(i) WPC does not accommodate a cutoff independent T-matrix at 0(Q 2 ) or 0(Q 3 ) when sub- 
leading potentials are treated as perturbations on top of the LO. 

(ii) RG invariance can be achieved by the modified power counting ([8]), based on modified naive 
dimensional analysis. This suggests that — 1/r 2 , the LO long-range potential in the model 



of Ref. [12|| , is not crucial for NDA to be applicable, 
(iii) At C((5 3 ), Iter-CD, Pert-CD and Iter- WPC show that in a limited range of cutoffs these three 



approaches produce similar phase shifts for 3 Pq, a conclusion similar to that of Refs. (itI . 



3. 



While it is instructive, we refrain from drawing the same conclusion for other channels; it 
may well be that the "common" window of cutoffs appears at different location for different 
channels. 



Ill ] has come to our attention. The conclu 



A simultaneous, coordinate-space calculation in Ref. 
sion drawn there for 3 Pq agrees with ours, that is, in agreement with NDA. However, Refs. 



10 



11] 



concluded a proliferation of six counterterms in each of the coupled channels, 3 Si — 3 Di and 
3 P2 — 3 F2, whereas NDA suggests three once WPC is corrected at LO. We defer to a further 



momentum-space calculation of the triplet channels 



261 ] that are subject to the singular attraction 



id, 



ii 



on the coupled channels can 



of OPE to investigate whether NDA or the conclusion of Refs. 
be verified. 
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Appendix A: 

In distorted wave expansion, the unitarity of the S'-matrix no longer rigorously holds. Therefore, 
the normalization of the T-matrix we adopt in the paper, T = — e 4<5 sin 5 / m^k, is no longer exact 
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at subleading orders. Suppose that, with vanishing, the T-matrix and the phase shifts have 
the following expansion: T = T(°) + T( 2 ) + T( 3 ) + • • • , 5 = ff(°) + ffW + 6® + .... Treating T^' ^ 
and ^( 2 )'( 3 ) as perturbations, one finds 

T (2), (3) = _ § (2), (3) e ^ 
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